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Abstract 

Molecular dynamics simulations are used to study the spreading of binary polymer nanodroplets in 
a cylindrical geometry. The polymers, described by the bead-spring model, spread on a flat surface 
with a surface-coupled Langevin thermostat to mimic the effects of a corrugated surface. Each droplet 
consists of chains of length 10 or 100 monomers with ~ 350 000 monomers total. The qualitative features 
of the spreading dynamics are presented for differences in chain length, surface interaction strength, 
and composition. When the components of the droplet differ only in surface interaction strength, the 
more strongly wetting component forms a monolayer film on the surface even when both materials are 
above or below the wetting transition. In the case where the only difference is the polymer chain length, 
the monolayer film beneath the droplet is composed of an equal amount of short chain and long chain 
monomers even when one component (the shorter chain length) is above the wetting transition and the 
other is not. The fraction of short and long chains in the precursor foot depends on whether or not both 
the short and long chains are in the wetting regime. Diluting the concentration of the strongly wetting 
component in a mixture with a weakly wetting component decreases the rate of diffusion of the wetting 
material to the surface and limits the spreading rate of the precursor foot, but the bulk spreading rate 
is not affected until the more strongly wetting component is removed completely. 

1 Introduction 

The spreading of liquid droplets on a surface is an important issue for several industries including adhesion, 
lubrication, coatings, and printing. Emerging nanotechnology in areas such as lithography and microfiuidics 
make the issue of droplet spreading on small length scales even more relevant. 

Just as blending bulk polymers can improve the physical properties of the resulting material, adding a 
second component to a spreading droplet can produce a desired change in surface tension or wettability. For 
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example, adding a surfactant to a droplet that ordinarily does not wet a surface can give a product that 
does wet the surface.^ In general, adding a second component provides more parameters with which to tune 
the material properties, but it also introduces more complex phenomena such as the interdiffusion of the two 
components and possible preferential wetting of one at either the liquid/solid or liquid/ vapor interface due 
to the difference in surface tension. 

Compared to homogeneous systems, the published literature on binary droplets is limited. Experiments 
on binary droplets have been concerned primarily with a nonwetting polymer solute in a wetting solvent. 
These experiments explored the concentration dependence of the equihbrium contact angle^ and found a 
leak-out transition^'"' where a film of pure solvent is in equilibrium with the droplet. Other experimental^ 
and theoretical papers^'^ focus on the equihbrium behavior of binary mixtures. The surface segregation of 
two-component polymer films has been studied using molecular dynamics (MD) simulation*"^" as well as 
numerical techniques based on Ising-Hke models. ^^"^^ Some MD simulations of droplet spreading have been 
performed for mixtures of solvent and oligomers^'''^^ with 4 000 monomers or mixtures of two oligomers^^ 
with up to 25 000 monomers, though even larger system sizes are required to adequately model the spreading 
dynamics. One work by the current authors explored reactive wetting of binary droplets in metallic systems. ^''^ 
Results therein demonstrated a difference in substrate dissolution rate with varying droplet composition, 
which in turn resulted in different wetting kinetics and equihbria. However, the combined phenomena of 
substrate dissolution with wetting prevents directly connecting those results to non-reactive wetting of binary 
polymer droplets. 

This paper presents molecular dynamics simulations of binary polymer nanodroplets differing in either 
polymer chain length or surface interaction strength. We explore the qualitative differences in the spreading 
behavior, surface monomer composition, and precursor foot composition for several systems. We also vary 
the relative concentration of the two components and analyze the velocity distributions of the spreading 
droplets. All of the simulations presented here are for droplets of ~ 350 000 monomers, which for single 
component droplets is sufficiently large to overcome finite size effects. ^^'^^ 

This paper is organized as follows. Section [2 presents a brief summary of the molecular dynamics 
simulation techniques and the application of the surface-coupled Langevin thermostat. It also describes 
the methods used to analyze the simulation results. The methods of analysis are described in more detail 
elsewhere. ^^'^^ SectionElcontains the simulation results for binary polymer nanodroplets for different polymer 
chain lengths and surface interaction strengths, and discusses the effects of varying the mole fraction of 
strongly wetting and weakly wetting components. Conclusions drawn from these simulations are presented 
in Sec. 21 
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2 Simulation Details 



Molecular dynamics simulations are performed using the standard coarse-grained model for polymer chains 
in which the polymer is represented by spherical beads of mass m attached by springs. We use a truncated 
Lennard-Jones (LJ) potential to describe the interaction between the monomers. The LJ potential is given 

by 
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where e and a are the LJ units of energy and length and the cutoff is set to r,; = 2.5 a. The monomer- 
monomer interaction e is used as the reference and all monomers have the same diameter a. For bonded 
monomers, we apply an additional potential where each bond is described by the finite extensible nonlinear 
elastic (FENE) potential^" with A; = 30 £ and Rq = 1.5 a. 

The substrate is modeled as a fiat surface since it was found previously^^ that with the proper choice of 
thermostat, the simulations using a fiat surface exhibit the same behavior as those using a realistic atomic 
substrate. Since simulating a realistic substrate requires several times the total number of atoms in the 
simulation, using the fiat surface greatly improves the computational efficiency. The interactions between 
the surface and the monomers in the droplet at a distance z from the surface are modeled using an integrated 
LJ potential. 
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where Syj is the monomer-surface interaction strength and Zc = 2.2a. 

We apply the Langevin thermostat to provide a reahstic representation of the transfer of energy in the 
droplet. The Langevin thermostat simulates a heat bath by adding Gaussian white noise and friction terms 
to the equation of motion. 



miTi = -AUi - maLri + Wi{t), (3) 

where rrij is the mass of monomer i, is the friction parameter for the Langevin thermostat, —AUi is 
the force acting on monomer i due to the potentials defined above, and Wi(t) is a Gaussian white noise 
term.^^ Coupling all of the monomers to the Langevin thermostat has the unphysical eflFect of screening the 
hydrodynamic interactions in the droplet and not damping the monomers near the surface stronger than 
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those in the bulk. To overcome this, we use a Langevin couphng term with a damping rate that decreases 
exponentially away from the substrate. We choose the form Jl{z) = 7£ exp {a — z) where 7£ is the surface 
Langevin coupling and z is the distance from the substrate. We generally use values of 7£ = 10.0 and 
3.0 for Eu, = 2.0 £ and 3.0 £, respectively, based on earlier work^^ matching the diffusion constant of the 
precursor foot for flat and atomistic substrates. The larger 7£ corresponds to an atomistic substrate with 
larger corrugation and hence larger dissipation and slower diffusion near the substrate. 

All of the droplets presented here are modeled as hemicylinders as described previously.^^ The droplets 
spread in the x direction and each system is periodic in the y direction with length Ly = 40 a and open 
in the other two directions. This allows a larger droplet radius to be studied in the cyhndrical geometry 
than in the spherical geometry using the same number of monomers. The binary droplets consist of either 
mixtures of chain length A'' = 10 and N = 100 polymers with the same surface interaction strength Sw or 
mixtures of chain length N = 10 with varying Ew Unless explicitly stated, the binary droplets contain an 
equal number of monomers of each component with initial droplet radius of Ro = 80 a for a total size of 
350 000 monomers. These droplets are large enough that finite droplet size effects are minimal.^* 

The equations of motion are integrated using a velocity- Verlet algorithm. We use a time step of At = 

1 /2 

0.01 T where t = a . The simulations are performed at a temperature T = e/ks using the lammps 

code.^^ Most of the simulations were run on 48 to 64 processors of Sandia's ICC Intel Xeon cluster. One 
milhon steps for a wetting drop of 350000 monomers takes 37 hours on 48 processors. 

For the simulations presented here, the instantaneous contact radii for the precursor foot and bulk regions 
are extracted every 400 r. The contact radius is calculated by defining a one-dimensional radial distribution 
function, as described previously.^* The precursor foot radius calculation includes all monomers that are 
within 1.5 cr of the surface whereas the bulk radius calculation uses monomers between 4.5 and 6.0 cr from 
the surface. 

3 Results and Discussion 

The additional degrees of freedom for binary droplets leads to a number of qualitatively distinct spreading 
characteristics. In the following figures, we show snapshots of several different droplets as they wet the 
substrate. In each case, the droplets start in roughly the same configuration with a contact angle near 
90°, but they each show differences in the precursor foot composition and spreading rate as well as the 
composition of the first layer above the substrate. 
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Figure 1: Profiles for binary droplets composed of equal mixtures of chain length = 10 with — 1-0 e 
(yellow) and e„ = 2.0e (blue) at t = 40 000 t (top) and t = 80 000 r (middle). Profile for = 2.0e (yellow) 
and Ew = 3.0 e (blue) at t = 40 000 r (bottom). The substrate (red) has a length of 600 tr and a width of 
40 cr in each profile. 7£ = 10.0 r~^. 

3.1 Substrate Interaction Strength 

The top two droplets in Fig. ^show a profile of a spreading droplet for the mixture of chain length A?^ = 10 in 
which half of the chains interact with the substrate with — 1.0 e and half with — 2.0 e for two different 
times. Previously we have shown that the wetting transition for droplets of chain length = 10 occurs at 

~ 1.75 e. This droplet is a mixture of a wetting and nonwetting polymer.^* For N = 10, = 1.0 e has 
a finite contact angle of 6*0 = 90°. At early times, a monolayer of the wetting component forms at the solid 
interface and wets the substrate. However, unlike the case of a homogeneous droplet in which all monomers 
interact with the substrate with e^, = 2.0 e,^^ the precursor foot in this case does not separate from the main 
droplet. Both components of the droplet follow the monolayer as it wets the substrate, but the spreading 
rate of the droplet is limited by the spreading rate of the monolayer. The evolution of the contact angle 
for the binary droplet with — 1.0 £ and = 2.0 £ is compared to the single component droplet with 
£u, = 2.0 £ in Fig. El In this case, adding the lower component slightly decreases the spreading rate of 
the droplet. This decrease is sufficient to produce a better fit to the hydrodynamic spreading model^® as 
compared to the pure Sw — 2.0 £ system. 

Different behavior is observed when both components wet the surface as shown in the bottom panel of 
Fig. n]for a mixture with £„, = 2.0 £ and 3.0 £. Despite the fact that both components wet the substrate, 
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Figure 2: Dynamic contact angle for the binary systems with an equal mixture of £w = 0.5 e and = 1-5 e 
(triangles) and an equal mixture of = 1-0 e and £w = 2.0 £ (squares) compared to the single component 
droplets with e^j = 1-5 e (diamonds) and = 2.0 £ (circles). N = 10, 7£ = 10.0 r""'^. 

a monolayer of the Sw = 3.0 £ component still forms at the solid interface since that is the strongly favored 
interaction. Although earlier simulations showed an accumulation of the more strongly wetting component 
in the first layer above the substrate/® here the first layer consists entirely of the more strongly wetting 
component. Unlike the = 1.0£/2.0£ case, the spreading rate of the precursor foot is fast enough to allow 
it to advance well ahead of the bulk of the droplet similar to the leak-out behavior observed experimentally 
for polymer/solvent mixtures.^''' A depletion region forms near the edges of the bulk droplet because the 
diffusion rate of the wetting component in the bulk is insufficient to replace the material forming the precursor 
foot. In Section we return to this case for varying composition of the two components. 

We also consider the case of a droplet containing a mixture of two nonwetting components, Ew — 0.5 £ 
and £u, = 1.5£ for iV = 10. As in the previous cases, a monolayer of the more wetting component forms 
between the droplet and the substrate. In this case, the droplet reaches an equilibrium contact angle of 
6*0 ^ 64° after 20 000 r. The evolution of the contact angle for this nonwetting system is compared to the 
Ew — 1.5 £ single component system in Figure IS The single component droplets have equilibrium contact 
angles of 6*0 — 120° for Sw = 0.5 £ and 44° for Sw = 1.5 £. Thus the contact angle of the binary system is 
more strongly infiuenced by the material that forms a monolayer on the substrate. 

3.2 Polymer Chain Length 

For droplets composed of two chain lengths of the same type of polymer, the spreading behavior is qualita- 
tively different. The top two droplets in Fig. El show profiles for the mixture of chain lengths = 10 and 
100 with — 2.0 £ and 7£ = 10.0 r^^. Although £^ = 2.0 £ for both components, the longer chains are in 



6 



Figure 3: Profiles for binary droplets composed of equal number of monomers of chain length = 10 (blue) 
and N = 100 (yellow) polymers with = 2.0 e and Yl = 10.0 t^^ at t = 40 000 r (top) and t = 80 000 r 
(middle). Same as above with Ew = 3.0 e and 7£ = 3.0 taken at t = 40 000 t (bottom). The substrate 
(red) has a length of 800 a and a width of 40 a in each profile. 

the nonwetting regime (e^ ~ 2.25 e)^^ whereas the shorter chains are slightly above the wetting transition 
(e^ ~ 1.75 e). As a result, a slow precursor foot consisting of approximately 80% iV = 10 monomers and 
20% N — 100 monomers extends from the bulk region. This is somewhat similar to the "leak-out" behavior 
observed experimentally in mixtures of polymer and solvent where the precursor foot consists entirely of 
solvent.^''' Unlike Fig. ^ a monolayer of the wetting component does not form on the substrate below 
the bulk region of the droplet. Instead, the first monolayer at the substrate beneath the droplet consists 
of an equal fraction of the two components. This can be understood by noting that the longer chains are 
nonwetting in this case due to their greater surface tension. The chains that are buried beneath the droplet 
feel no infiuence of the liquid/vapor interface, so there is no preference as to which chain length is in con- 
tact with the substrate. The chains at the upper surface of the droplet, both in the bulk region and the 
precursor foot, are infiuenced by the liquid/ vapor interface and there one finds an abundance of the shorter 
chains. The liquid/vapor surface tensions for the = 10 and N — 100 chains were previously found to 
be 0.84 and 0.96 ± 0.02 e/a^, respectively.^^ The surface tension for the 10/100 mixture was found to be 
0.90±0.02 e/(T^ indicating a roughly equal infiuence by each component on the surface tension of the mixture. 
Thus, the surface tension calculation is not precise enough to account for the abundance of shorter chains 
at the liquid/vapor interface. 

If Ew is increased so that both chain lengths are in the wetting regime, the behavior is similar to that of 
a single-component droplet. This is shown in the bottom of Fig. Elfor the mixture of chain lengths = 10 
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Figure 4: Velocity profiles of half of the droplet showing the iV = 10 chain length (a), and the N — 100 
chain length (b), for the binary droplet with Sw = 3.0 e and 7£ — 3.0 at t = 20 000 r. 

and 100 where Sw = 3.0 e and 7£ = 3.0 r~^. Here, both chain lengths are in the wetting regime and the 
rapidly spreading precursor foot is composed of an equal number of monomers of the two chain lengths. One 
consequence of this is the N — 100 chains diffuse across the substrate more rapidly when mixed with TV = 10 
chains than in the pure N — 100 droplet. Although the longer chains have a higher surface tension, the 
equilibrium state for both components is a 0° contact angle so there is no segregation of the two components 
in the precursor foot as this droplet spreads. The lack of segregation for a mixture of two chain lengths with 
the same Sw was previously reported by Voue et al. for chain lengths 8 and 16.^^ 

The droplet velocity profiles provides another method to analyze the dynamics of the spreading droplet. 
Figure 01 shows the velocity profiles for the mixture of chain lengths 10 and 100. Although both components 
show the same spreading pathway, the shorter chain lengths tend to move slightly faster at the liquid / vapor 
and liquid/solid interfaces. 

3.3 Varying Composition 

The ability to form a film from a non- or weakly wetting polymer by adding a wetting component is of 
significant practical interest. To analyze the effects of the concentration of the two components on the 
spreading dynamics, we simulate droplets containing a strongly wetting component {Ew — 3.0 e) at different 
mole fractions, x^et — 0.10, 0.25, and 0.50 with a weakly wetting component {Ew = 2.0 e). Since = 10 
for both, the only difference in the two components is the surface interaction. Profiles of the spreading 
droplets after 60 000 r are shown in Figure |3 In each case, a monolayer of the strongly wetting component 
forms on the surface and spreads outward. The rest of the drop then spreads on this layer. However, for 
Xwet = 0.10, the surface monolayer never becomes fully homogeneous as a segment of the weakly wetting 
component remains at the edge of the foot and continues to get pushed outward as the droplet spreads. The 
size of the depletion region at the droplet edge increases as Xwet decreases since less of the strongly wetting 



(a) 



w j K. H i'. i'/t l\ V\ \ 



V N vp. v^^^^^^^^^|^ i t \ ti t t t t 



(b) 



8 



Figure 5: Profiles of three binary droplets of chain length N = 10 with = 2.0 £ (yellow) and Sw = 3.0 e 
(blue). The mole fractions of the strongly wetting component are x^et = 0.10 (top), 0.25 (middle), and 0.50 
(bottom), t = 60 000 r, 7£ = 10.0 r~^. The substrate (red) has a length of 800 a and a width of 40 a in each 
profile. 

polymer is available to form the precursor foot. The rate of spreading of the precursor foot decreases as 
Xwet decreases because the supply rate of the strongly wetting material from the bulk to the surface limits 
the spreading rate of the precursor foot. This is more clearly demonstrated in Fig. El where the contact 
radii of both the foot and bulk regions are shown as a function of time for each of the three compositions 
and compared to previous results^® for Xwet = 0. Figure shows that the precursor foot spreading rate 
decreases as Xwet is reduced. The effective diffusion constant of the precursor foot, D^ff — 5-^, decreases 
from 4.7 to 3.8, 2.6 and 0.6it^/t as Xwet decreases from 0.50 to 0.25, 0.10 and 0.0, respectively. Figure 

shows no significant effect of x^et on the bulk droplet spreading rate for the duration of the simulation 
for Xwet — 0.5 to x^et =0.1. However, there is a strong effect on the spreading rate when increasing the 
concentration from x^et — 0.0 to x^^et =0.1. In all of these cases, the bulk contact radius spreads as with 
X = 1/5,^^ which is expected from the kinetic model for spreading in the cylindrical geometry used here. 

The velocity profiles of the strongly wetting and weakly wetting components at t = 70 000 r are shown 
in Fig. dfor Xwet = 0.25. After the foot has extended and a depletion region has formed at the edge of 
the droplet. Fig. shows that the weakly wetting material moves down from the liquid/vapor surface 
and up from the solid/liquid surface to allow the monolayer of the strongly wetting material to form. Very 
little motion is seen in the center of the droplet. In Fig. [3d, the strongly wetting material shows similar 
behavior at the liquid/vapor surface, but the data is noisier due to the smaller mole fraction of strongly 
wetting material. This material moves through the depletion region at the edge of the droplet and onto 



9 



6.0 
*o 4.0 



(^"2.0 



0.0 
1.5 

o 



OS 



»0.5 



(a) 



(b) 



0„^-^ 



_i_ 



20000 



40000 
time [x] 



60000 



80000 



Figure 6: Contact radii of the precursor foot (a) and bulk (b) for the droplets composed of x^et = 
(dotted curves), x^et = 0.10 (dashed curves), 0.25 (dot-dashed curves), and 0.50 (dot-dash-dashed curves). 
Sw = 2.0 e for the weakly wetting component and = 3.0 e for the strongly wetting component. N = 10, 
Yl = 10.0 r-^ 

the precursor foot where it rapidly diffuses outward. The same qualitative behavior is observed for different 
values of Xwet, though the data becomes increasingly noisy as x^et decreases. 



4 Conclusions 

The spreading dynamics of binary polymer nanodroplets are studied using molecular dynamics simulation. 
We demonstrate qualitative differences in the spreading behavior of binary droplets due to differences in 
surface interaction strength, polymer chain length, and composition. 

When the two droplet components differ in surface interaction strength, the more strongly wetting com- 
ponent forms a monolayer film on the surface even when both materials are either above or below the wetting 
transition. For cases that also produce a rapidly spreading precursor foot, a depletion region of the more 
strongly wetting component forms starting at the edge of the droplet since the interdiffusion rate in the 
droplet is slower than the rate at which material is withdrawn into the precursor foot. 

For differences in the polymer chain length, the monolayer film beneath the droplet is composed of an 
equal amount of short chain and long chain monomers. This is true even when one component (the shorter 
chain length) is above the wetting transition and the other is not. In this case, the precursor foot is composed 
primarily of the wetting component. When both components are above the wetting transition, the surface 
monolayer and precursor foot are both composed of an equal amount of short and long chain monomers. 

The formation of a monolayer of the more strongly wetting component is studied in greater detail by 
considering different compositions of the mixture of strongly wetting and weakly wetting components. For 
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Figure 7: Velocity profiles of half of the droplet showing the weakly wetting component, Syj = 2.0 e (a), and 
the strongly wetting component, = 3.0 e (b), for the x^et = 0.25 binary droplet of chain length A'' = 10 
with Yl = 10.0 at t = 70 000 r. 

each case, a monolayer of the strongly wetting component forms on the surface and spreads outward. The 
spreading rate of this precursor foot decreases as the concentration of the strongly wetting component 
is reduced, but the spreading rate of the bulk droplet shows no significant dependence on the droplet 
composition until all of the more strongly wetting component is removed. Velocity profiles show that the 
material for the precursor foot is supplied by the wetting material near the hquid/ vapor surface as it diffuses 
through the depletion region and onto the substrate. 
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